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We measure the particle-level forward-backward production asymmetry in bb pairs with masses 
(m^j) larger than 150 GeV/c^, using events with hadronic jets and employing jet charge to distinguish 
b from b. The measurement uses 9.5 fb~^ of pp collisions at a center of mass energy of 1.96 TeV 
recorded by the CDF II detector. The asymmetry as a function of is consistent with zero, as well 
as with the predictions of the standard model. The measurement disfavors a simple model including 
an axigluon with a mass of 200 GeV fc? whereas a model containing a heavier 345 GeV/axigluon 
is not excluded. 

PACS numbers: 14.65.Ha 


I. INTRODUCTION 

In recent years, the values of the forward-backward 
asymmetry (Ape) of top-quark-pair production measured 
at the Tevatron proton-antiproton collider at Fermilab 
have been consistently larger [1] than those predicted by 
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TABLE I. Predicted values of Ape in the standard model [11] 
and in two models with an axigluon [9]. In the second column, 
the first contribution to the uncertainty is due to neglected 
higher-order terms, and the second contribution is due to 
the combined effect of varying the factorization and renor¬ 
malization scales. The assumed axigluon mass is listed, and 
its assumed width is 25 % of the mass in both cases. The 
selection requirements imposed on all three calculations match 
the event selection requirements employed in our analysis: 
pseudorapidity <1.1. 


mt,i range 

A^Bibb) [%] 


[GeV/c^] 

SM 

200 GeV/c^ 

345 GeV/c^ 

[150,225] 

2.43±0.73to.oi 


_1 Q + 0'4 
-*-•^-0.9 

[225,325] 

4.61T 1.38t|];]^ 

20.4l?;o 

_q q+O-6 

[325,1960] 8.70 ± 2 . 61 ^ 0 . 5 ] 

20.2tli 

16.41°;^ 


the standard model (SM) [2-7]. Further study of this 
phenomenon has led to a number of proposed extensions 
of the SM [8]. One specific class of such models is the 
low-mass axigluon [9]. These models include a new axial- 
vector boson with a mass below the tt threshold and a 
natural width broad enough to evade detection in light- 
quark resonance searches, and predict a non-zero forward- 
backward asymmetry due to the interference between 
amplitues mediated by the gluon and the axigluon. In 
the simplest models, the axigluon has equal couplings to 
all the quarks. Hence, a stringent test of such models is 
the measurement of the forward-backward asymmetry of 
pair production of other quark flavors, such as bottom. 

The forward-backward asymmetry for fermion- 
antifermion production is defined as 


Afb = 


np — riB 

np + ub' 


( 1 ) 


where np is the number of events where the fermion is 
forward of the antifermion in rapidity [10] {Ay = yb —Hi > 
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0) and nB the number where it is backward (Ay < 0). 
This definition is invariant under boosts along the beam 
axis. 

At hadron colliders, bb pairs are almost exclusively 
produced by the strong interaction of quarks and gluons 
(QCD). The vast majority of bb pairs are produced via 
gluon-gluon fusion, which yields no asymmetry due to 
the symmetric initial state. In the qq —> bb process, a 
positive ApBibb) arises from higher-order QCD corrections 
involving either real or virtual gluons. Since the valence- 
quark parton-density functions dominate over the gluon 
parton-density functions at large Bjorken x, we enhance 
the contribution of the quark-antiquark initial state over 
the symmetric gluon-fusion background by requiring large 
values of bb mass, mf,j, > 150GeV/c^. 

A number of theoretical predictions of the SM value 
for ApB{bb) at high mass have been reported, using var¬ 
ious techniques and kinematic requirements and yield¬ 
ing a range of predictions [4, 7, 11]. Kiihn and Ro¬ 
drigo [4] computed that the asymmetry for bb pairs with 
> 300 GeV and production angle |cos0*| < 0.9 falls 
in the range 4.3_% to 5.1%. Manohar and Trott [7] 
computed Ape (65) in several bins of m^f,. They found 
Apnibb) = 0.4% inclusively and Apsibb) = 7.8% to 8.1 % 
for 350 < < 650GeV/c^, depending on the choice of 

the factorization and renormalization scales. Grinstein 
and Murphy [11] compute predictions using a variety 
of kinematic requirements, including an estimate using 
our kinematic requirements. The results of this calcula¬ 
tion [11] are summarized in Table I. 

Models containing a low-mass axigluon that has the 
same couplings to all quark flavours also predict a forward- 
backward asymmetry in bottom-quark-pair production. 
This asymmetry arises from the interference between the 
gluon and axigluon in the same fashion as the Z in¬ 
terference of the Drell-Yan process [12]. One key feature 
of this interference is that the App changes sign at the 
mass pole of the heavy axial resonance [9]. Observation 
of such a sign flip would provide significant indication of 
the underlying dynamics. The predictions of two different 
representative axigluon models are given in Table I [9]. A 
relatively light 200 GeV /(? axigluon has been previously 
studied in the context of the top-quark asymmetry [13], 
and a heavier 345 GeV /(? axigluon is of interest because 
its mass is just below the tt threshold. Both of these 
axigluons are assumed to have a width equal to 25 % of 
their mass. 

The predictions of both the SM and the models of 
non-SM physics under consideration are estimated at the 
parton level, and do not explicitly include the effects of 
hadronization. 

Although the LHCb experiment has measured a re¬ 
lated quantity [14], App can only be measured in proton- 
antiproton collisions [15]. At the Tevatron, the DO ex¬ 
periment has measured the forward-backward asymmetry 
of meson production [16]. The momenta of the 
mesons in that measurement constrain the masses of the 
bb pairs to be smaller than the range probed in this mea¬ 


surement. 


II. EXPERIMENTAL APPARATUS AND 
EVENT SELECTION 


The Fermilab Tevatron is a proton-antiproton collider 
with a center-of-mass energy ^/s = 1.96 TeV. The CDF II 
detector [17] is an azimuthally symmetric magnetic spec¬ 
trometer with a large tracking volume inside a solenoid. 
Outside the solenoid are sampling calorimeters. The 
calorimeters are further surrounded by the steel flux re¬ 
turn of the solenoid and muon detectors [17]. 

The data used for this analysis were collected by the 
CDF II detector using three online event selections (trig¬ 
gers), which require at least one jet with transverse energy 
Et > 50, 70, or 100 GeV, respectively [10]. To control the 
trigger rate, the two lower-threshold triggers only accept 
one event out of every 100 or 8 events that satisfy their 
requirements, respectively, while the highest-threshold 
trigger accepts all events meeting its requirements. After 
data-quality requirements and trigger acceptance rates, 
the integrated luminosities of the samples are 95pb“^, 
1.2fb“^, and 9.5fb“^, respectively. 

The offline event selection requires at least two jets [18, 
19] with Et > 20 GeV and pseudorapidity \r]\ < 1.1. 
Of the jets that pass these requirements, exactly two 
must contain a secondary vertex (separated transversely 
from the primary pp interaction vertex) consistent with 
the decay of a 5 quark (6-tagged jets), identified using 
a secondary-vertex identification algorithm [20]. The 
observed invariant mass of the two 6-tagged jets is 
required to be at least 150GeV/c^. 

In order to maximize the statistical significance of the 
analysis, we separate the data into several subsamples 
depending on (three subsamples, see Sec. Ill), on the 
estimated charges of the partons that produced the jets 
(four subsamples, see Sec. Ill A), and on the quality of 
the b tags (three subsamples, see Sec. IIIC), yielding a 
total of 36 subsamples. 


III. METHODOLOGY 


We discriminate jets originated from b quarks {b jets) 
from jets originated from b quarks (6 jets) using the 
momentum-weighted average of the charges of the parti¬ 
cles associated with each jet (jet charge), 




( 2 ) 


where the sum i is over all tracks in the jet, qi is the 
charge of the corresponding particle, and pi and pjet are 
the momentum vectors of the particle and of the jet, 
respectively. The exponent 0.5 was chosen to maximize 
the power of the jet charge in separating b jets from b 
jets [21]. 
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The asymmetry depends on the mass of the bb pair, 
both in the SM and in models with an axigluon. To study 
the behavior of Ape as a function of we divide the 
sample into several ranges (bins) of The choice of 
bins is motivated by the trigger. Each jet-energy-trigger 
threshold efficiently selects events only over a limited 
range of Events with 150 < < 225GeV/c^ are 

selected with a jet transverse-energy threshold of 50 GeV, 
and the thresholds of 70 and 100 GeV are used to select 
events with 771^5 in the ranges 225 < < 325 GeV jt? 

and 325GeV/c^ < respectively. 

The measurement of A-p-aibb) must account for non- 
bb backgrounds, detector effects, and the dilution of the 
asymmetry due to misidentification of b jets as b jets and 
vice versa. We therefore define a detector-level Ape, 

det ^ 1 (iVp - - {Nb - 

2P - 1 iVp - ATbkgd +Nb- ’ 

where Np and Nb are the observed numbers of events, 
and are the estimated numbers of back¬ 

ground events, and P is the probability to make the 
charge assignment correctly. To account for additional 
effects, such as the finite resolution and acceptance of 
the GDF II detector, we employ a Bayesian technique to 
measure the Tpp at particle level. “Particle level” refers 
to quantities reconstructed from final-state, color-confined 
particles with lifetimes greater than 10 ps [ 22 ]. 

The background levels are determined with a data- 
driven technique based on the ^-tagged vertex mass as 
described in Sec. IIIC. The calibration of the charge 
misidentification and the calculation of P is presented in 
Sec. Ill A, and the background asymmetries are discussed 
in Sec. HID. Finally, we use a Markov-chain Monte Carlo 
calculation to derive particle-level results by identifying 
the maxima of the marginalized posterior probability 
densities and by constructing associated credible inter¬ 
vals. The correction to the particle level, including the 
effects of mismeasurement and acceptance, are described 
in Sec. IIIG. 


A. Identification of the b jet using jet charge 

The forward or backward assignment is performed using 
the momentum-weighted track charge, or “jet charge” (see 
Eq. (2)), for each of the two 6 -tagged jets. Distributions of 
the jet charge (Qjet) for 6 jets and & jets are shown in Fig. 1. 
We use the difference of the two jet charges AQ = Q 1 — Q 2 
to make the assignment: if AQ is negative, the jet with 
charge Qi is considered to be the b jet (because the b 
quark is negatively charged), and if AQ is positive, the 
jet with charge Q 2 is considered to be the b jet. 

The performance of the jet-charge algorithm is cali¬ 
brated with data, by dividing the distribution of Qjet into 
bins. The bin edges are —0.25, 0, and 0.25, for a total of 
four bins (the lowest- and highest-charge bins are open). 
We arbitrarily assign Qjet values of —0.5, —0.25, 0.25, 



FIG. 1. Distributions of the charge of b and b jets in Monte 
Carlo simulation of two-jet events. 


and 0.5 to the jets falling into each respective bin. This 
translates into five bins of |AQ| with values 0, 0.25, 0.5, 
0.75, and I.O. The bin with |AQ| = 0 is not informative 
because it gives no indication of how to assign the 6 jet, 
so only four bins of |AQ| are used. 


B. Charge-identification probability 


The probability P of correctly assigning the b jet and 
b jet, introduced in Eq. (3), is estimated from the data. 
We evaluate this in each bin of |AQ|, so that 


7'o.5(1—^0.25) + (1 ~7'0.5)P0.25’ 
Po.25 


Ab.25 = 

^0.5 = 

-P 0.75 = 

Pl.o = 


^0.5(1 — ^0.25) 


Po.25 + (1 -Po. 25 )^’ 

P 0 . 25 P 0.5 


P 0 . 25 P 0.5 + (1 — P0.25)(l ~ P 0 . 5 ) ’ 

PI 5 


and 


PI5 + (1 -Po.5)^’ 


(4) 


where the subscript of P indicates |AQ|, and the various 
P\Qiet\ the probabilities that a jet with a binned charge 
of I Qjet I is correctly identified as 6 or & by the sign of Qjet- 
These expressions are derived by exhaustively considering 
all the possibilities. For example, |AQ| = 1.0 implies that 
one jet has Qjet = -1-0.5 and the other jet has Qjet = —0.5. 
If the b jet has charge —0.5, then the use of AQ allows for 
a correct assignment of the b and b jets. This case occurs 
with probability Pq 5 , because both jets (with jQjetl = 0.5) 
have a jet charge whose sign matches the sign of the charge 
of the originating quark. The other possibility is that 
we misidentify both jets, which occurs with probability 
(1 — Po.s)^- The denominator in Eq. (4) is the sum of 
these two probabilities, and the numerator only contains 
the correct case. 

In order to measure P|Qjj,th measure the number of 
opposite-charge events, AqCi number of same-charge 
events, AsC; in the data in each subsample. These 
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numbers are corrected for the presence of background 
using the calibrated 6-fractions (see Sec. Ill C) and the 
background model described in Sec. HID. We assume 
that the remaining sample is composed of 66 with no 
contamination from 66 or 66, and compute the opposite- 
charge fraction Foe = Noc/{Noc + iVsc)- 

In events in which both jets have |Qjet| = 0.25, the 
opposite-sign fraction is expressed as 


between high-momentum particles and the detector ma¬ 
terial. This results in an excess of positive over negative 
false tags. We exploit this relationship between the num¬ 
ber of positively- and negatively-tagged light-flavor jets 
and compute the number of true 66 events using 

N,-, = ^(iV++ - \N+_ + (8) 


^’Sc“=PL5 + (l-P0.25f, (5) 

where the term 25 arises from events in which both 6 
jets have the correct sign, and the term (1 — po. 25 )^ arises 
from events in which both 6 jets have the wrong sign. 
Similarly, we have 

pO.25 ° ® = pq 25 P 0.5 + (1 ~ Po. 25 )(l ~ P 0 . 5 ): and (6) 
F°-t°-^=pl, + {l-po.,r, (7) 

for events containing one jet with IQjetl = 0.25 and one 
jet with IQjetl = 0.5, and for events in which both jets 
have IQjetl = 0.5, respectively. We measure each Fqc 
value in the data and solve for each and then for 

each of the four P\aq\- 

C. Sample purity 

Equation (3) also requires knowledge of the rate of 
background events, and We obtain the 

background yields by estimating the number of 66 events 
and subtracting it from the number of events ob¬ 
served. 

We divide the data into subsamples with varying 6-tag 
quality, which in turn provides subsamples of varying 
66 purity and improves the statistical power of the mea¬ 
surement. Since we apply 6-tags to jets based on the 
presence of a secondary vertex, the quality of the 6-tag 
is based on the confidence level of the identification of 
the secondary vertex. Specifically, the significance is 
based on the distance L 2 D, in the plane perpendicular 
to the beam, between the primary and secondary ver¬ 
tices, projected onto the jet momentum [20]. Jets with 
a significance (|F 2 d| /o'(7^2d)) greater than 20 are high 
significance or “H” tags, and jets with a lower significance 
are “L” tags. This results in the following independent 
subsamples ranked in order of increasing 66 purity: LL 
(both jets are L tagged), LH (one L and one H tag), and 
HH (both H tags). 

We estimate the number of events in which both 6- 
tagged jets are genuine 6 jets by counting events in which 
one or both of the 6-tagged jets have a negative F 2 D- 
These “negative” tags are predominantly false tags from 
light-flavor jets and are a consequence of the finite po¬ 
sition resolution of the tracking system. We expect the 
rate of false tags from this source to be equal for positive 
and negative tags. There are additional false positive tag 
contributions from hyperon decays and from interactions 


where iV+_|_ is the number of observed positive double¬ 
tag events, _is the number of events with one of 

the tags negative, N _is the number with both tags 

negative, A is the ratio of positive to negative false tags, 
and ^ = 1 — 2Ar -|- A^r^ is a factor that accounts for the 
presence of negatively-tagged 6 jets. Here, r = e^/e)]' is 
the ratio of the negative to positive tag efficiencies for 6 
jets, estimated using Monte Carlo simulation and data. 
Finally, is corrected by a scale factor derived from 
MC to account for a bias resulting from the presence of 
charm jets [23]. 
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FIG. 2. Distribution of net mvTX, for H (a) and L (b) tagged 
data, with fit overlaid. The negative tag rates are subtracted 
from the positive tag rates. The histograms are stacked on 
top of one another. 


We measure the ratio of positive to negative light-flavor 
tag rates A following Ref. [23]. Since the MC indicates no 
strong dependence of A on jet Et, we use the same value 
for all three dijet mass bins. We study the distribution of 
the invariant mass of all charged particles associated with 
the secondary vertex, myxx [20], for both negative and 
positive tags. We subtract the negative tag distribution 
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FIG. 3. Estimated asymmetry of the backgrounds in each 
bin with tag type LL, as a function of |AQ|. 
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from the positive tag distribution, yielding the “net” dis¬ 
tribution. The relative rates of MC-derived templates for 
6 , c, and light-flavor jets, left as free parameters, are fit in 
the net tovtx distribution observed in data (see Fig. 2). 

These net flavor fractions allow for the normalization of 
the complete positive- and negative-tagged distributions 
of mvTX for all jet flavors. We find that the negative- 
tagged component of the MC templates must be scaled 
by 1.5 (1.3) to match the total number of low (high) 
significance negative tags observed in data. This scal¬ 
ing is propagated into the positive-tagged region of the 
templates as this population of “fake” tags is considered 
to be positive-negative symmetric, leaving the original 
“net” distributions unchanged [20]. From these corrected 
templates, we then obtain A for low- and high-significance 
6 -tagged jets. We also correct the values of r used in the 
three dijet mass bins for the higher negative tag rates 
observed in the data. 


FIG. 4. Smearing matrix relating measured dijet mass to 
particle-level dijet mass. The matrix is evaluated from MC, 
and depends on |AQ|, but it is shown summed over bins of 

|AQ|. 


asymmetry of the background in the signal sample (see 
Fig. 3). 

The asymmetry of the background is negative, consis¬ 
tent with u-u scattering as the most important component. 
The u quark tends to follow the incoming proton direction 
in t-channel scattering which would seem to yield a posi¬ 
tive ^fb, but the opposite charge of the u and b quarks, 
combined with the definition of the asymmetry in terms 
of the jet charge, reverses the sign of the asymmetry. 

Together with the background yields, the estimated 
background asymmetries are used to estimate the number 
of forward and backward background events, and 

Eq. (3). 


D. Background model 

The non-bb background includes a wide variety of 
physics processes, such as 6-|-mistag, uu, dd, gluon jets, 
etc. Of primary concern is the scattering of the valence 
quarks, which proceeds through the t-channel exchange of 
a gluon and exhibits a forward-backward asymmetry due 
to the forward Rutherford peak. Rather than attempting 
to tune the simulation to reproduce all of these processes 
in the proper amounts, we use data in a control region that 
is expected to be enhanced in background and depleted 
of bb signal. The control region is defined by loosening 
the requirements of the 6-tagging algorithm (looser than 
the low-significance “L” tag), and requiring that at least 
one of the tagged jets is negatively tagged. The sam¬ 
ple of events satisfying these criteria is non-overlapping 
with the signal sample, and is composed almost entirely 
of background but is kinematically very similar to the 
signal sample. The forward-backward asymmetry of this 
sample is assumed to approximate the forward-backward 


E. Jet energy mismeasurement 

Although we correct the energy of each jet for known 
effects [19], the energy of a measured jet does not ex¬ 
actly match the energy of the corresponding particle-level 
jet [22]. Mismeasurement of jet energies principally affects 
the measurement of the dijet mass. Since we use wide 
mass bins, this effect is small. We estimate the effect using 
MC samples of dijet events produced with the pythia 
event generator at leading order in the strong coupling 
constant, using the CTEQ5l parton distribution func¬ 
tions. We select events that have a pair of bottom quarks, 
whether produced directly or in the parton shower, to 
produce the signal model. 

This signal model allows the determination of a matrix 
relating the measured dijet mass to the particle-level dijet 
mass, as a function of |AQj. The matrix, summed over 
[AQ], is shown in Fig. 4. In 80 to 95% of events, the 
dijet mass is reconstructed in the same bin. The relative 
uncertainty on each element of the matrix ranges from 
a few percent on the diagonal up to very large values 
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TABLE II. Numbers of events observed in the various subsam¬ 
ples. 


nifji range AQ tag type N{Ay > 0) N{Ay < 0) 

[150,225] 0.25 

LL 

483 

472 


LH 

774 

841 


HH 

465 

469 

0.5 

LL 

734 
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LH 

1223 

1240 


HH 

682 

691 

0.75 

LL 

460 

474 


LH 

854 

882 


HH 

483 
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1.0 
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FIG. 5. Detector-level dijet mass spectra in each of the three 
irifji subsamples, showing the spectrum resulting from the 
non-66 component only as well as from the sum of the bb and 
non-66 components (“Modelilo-”). The agreement between 
data and simulation is good, suggesting that the particle-level 
dijet mass in MC simulation and the smearing matrix are also 
correct within the uncertainties. The non-6& component is not 
included in the smearing matrix. 


F. Selection effects 


far off the diagonal. The fraction of events in which the 
detector-level dijet mass migrates to a different bin than 
that of the particle-level dijet mass necessitates a correc¬ 
tion, and increases the uncertainty in the measurement 
of the particle-level Ape as a function of dijet mass (see 
Sec. IIIG). The correction relies on the MC description of 
the particle-level dijet mass to be approximately correct. 
Figure 5 shows that the detector-level dijet mass is cor¬ 
rectly described by the MC, which supports the reliability 
of the smearing matrix. 


We also estimate the effect of the potentially asymmet¬ 
ric acceptance of the detector and analysis selection. We 
apply the analysis event selection to the simulated events 
to estimate the fractions of forward (ep) and backward 
(cb) events that pass the analysis selection. We perform 
this estimation in each bin of particle-level dijet mass. 
Only the ratio of the forward to the backward acceptance 
affects the analysis, so we compute the ratio R = ep/cB- 
We also estimate the uncertainty on this ratio by varying 
the renormalization and factorization scales used in the 
signal model by a factor of two and the jet-energy scale 
within its uncertainty. 
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G. Extraction of the particle-level asymmetry 

To infer the asymmetry at the particle level, we con¬ 
struct a Bayesian model to describe the data. This model 
combines all effects discussed in the preceding sections, 
propagates the uncertainties, and allows the data to con¬ 
strain the uncertainties when possible. The parameters 
of the model and their assumed prior probability distri¬ 
butions are as follows: 

1 - /mqt is the bb fraction in each bin of detector-level 
dijet mass M, charge difference Q, and tag-quality 
T. The prior probability distribution for / is a 
normal distribution centered at the calibrated value, 
with a width equal to the residual uncertainty from 
the calibration. 

2. Pmq is the correct charge probability described in 
Eq. (4). The priors for po.s and po .25 are normal 
distributions with mean and uncertainty taken from 
the calibration. 

3. and rates of forward and back¬ 

ward events in the background-dominated sideband. 
From these, we calculate the background asymmetry, 
which is assumed to be consistent with the asym¬ 
metry of the background in the signal region. The 
prior for each of these is the gamma distribution. 

4. J is the shift in the jet energy scale. We coherently 
shift all of the jet energies in every MC event by 
J times the jet-energy uncertainty [19]. The prior 
is a normal distribution with a mean of zero and 
standard deviation of one. 

5. Sm'mq{J) is a matrix describing the contribution 
of events with particle-level dijet mass M' to the 
various bins of measured dijet mass and charge 
difference. The matrix is a function of J. The 
prior is taken from the rate and uncertainty in 
simulation, and the matrix is normalized so that 

I]m' Sm'mq{J) = 1- 


9. is the bb asymmetry after acceptance and se¬ 
lection effects. It is a function of Am' and Rm', 

.acc ^ Bm'(1 + Am') - (1 - Am') 

^ Bm'(1 + ^M') + (1 ~ ^M') 


We compute the rate 0 of forward and backward events 
expected in data, and compare this rate to that observed 
in data (see Table II) via a Poisson likelihood, given an 
observation of k events, [L{9\k) = /kl. The rates 
are 


^Forward 

^MQT 


flvIQT 


1 + A--%2Pmq-1) 


Sm'mq{J) 


M' 


+ (1 - Imqt) 


pbkgd 
^ MQT 

pbkgd . nbkgd 
^ MQT “b ^MQT. 


X CFmQT, 


(9) 


and 


qBackward 

^MQT 


I - A^%2Pmq - 1) , . . 

JMQT 2_^ -^-‘^M'MQW j 


M' 


+ (1 - Imqt) 


nbkgd 

^MQT 

pbkgd . pbkgd 
^ MQT “b ^MQT. 


X (JmQT- 


( 10 ) 


The prior probability densities described above, together 
with this likelihood, fully specify the posterior probability 
density for the parameters. To estimate this posterior den¬ 
sity, we employ Markov-chain Monte Carlo sampling [24]. 
This technique provides us with samples from the pos¬ 
terior probability distribution over the parameter space. 
We marginalize the nuisance parameters and obtain the 
posterior density for Am', the asymmetry in each bin 
of particle-level bb mass. The marginal distributions are 
shown in Fig. 6 , and the marginalized samples are pro¬ 
vided in the supplemental material [25]. 


IV. RESULTS AND CONCLUSIONS 


6 . (JMQT is the rate of events in each bin of detector- 
level bb mass, charge difference, and tag quality. 
This parameter has a uniform prior over the non¬ 
negative range. This parameter is necessary because 
the simulation does not accurately predict the over¬ 
all event rate as a function of mass, charge difference, 
and tag quality. 

7. Rm' is the ratio of the forward to the backward 
acceptance. The prior is a normal distribution 
with mean and width taken from the calibration 
described in Sec. IIIF. 

8 . Am' is the bb asymmetry in bins of particle-level 
mass. This is the parameter we wish to measure. 
We use a uniform prior from [—1,1]. 


To characterize the posterior and describe the mea¬ 
surement, we find the highest probability-density credible 
intervals at 68 % and 95 % credibility for Am' in each 
particle-level mass bin. The posterior densities, along 
with the intervals describing them, are shown in Fig. 6 . 
The red vertical bands (with solid lines) represent the 
theoretical predictions from the SM [II], while the blue 
and magenta bands (with dashed and dotted lines, respec¬ 
tively) represent the predictions from two axigluon models 
(see Table I). Although the measurement is performed at 
the particle level, the predictions are at the parton level 
and do not include the effects of hadronization. Because 
hadronization effects are not expected to be large, we do 
not hesitate to interpret the results. 

The measured asymmetries, summarized in Fig. 7, are 
—7.4lg;g%, and 6 .ll)^; 7 % in the low, middle. 
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and high mi,i, bins, respectively. These results, which 
account for the effects of backgrounds, charge misiden- 
tification, detector resolution, and nonuniform detector 
acceptance, are consistent with zero and with the stan¬ 
dard model prediction [11] in each bin. Only 0.24% of 
the posterior probability density in the middle mass bin 
has an Ape larger than predicted by the lighter axigluon 
model [9] with a mass of 200GeV/c^. Accounting for 
the look-elsewhere effect following Ref. [26], this is suffi¬ 
cient to exclude the lighter axigluon at more than 95 %. 
The measurement is unable to exclude the heavier ax¬ 
igluon with a mass of 345GeV/c^. This measurement 
reduces the allowed parameter space for light axigluon 
models used to explain the top-quark forward-backward 
asymmetry. 


FIG. 6. Marginal posterior probability distribution of asym¬ 
metry in each bin of particle-level bb mass. The inner and 
outer bands represent the 68 % and 95 % credible intervals, 
respectively. 
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FIG. 7. Values of maximum a posteriori signal asymmetry 
as a function of bb mass. The error bars represent the 68 % 
credible intervals. 
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